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Abstract 

The observation potential of the decay B + — * K + K + 7r~ with the 
ATLAS detector at LHC is described in this paper. In the Standard 
Model this decay mode is highly suppressed, while in models beyond 
the Standard Model it could be significantly enhanced. To improve the 
selection of the K + K + 7r~ final state, a charged hadron identification 
using Time-over- Threshold measurements in the ATLAS Transition 
Radiation Tracker was developed and used. 
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1 INTRODUCTION 



1 Introduction 

Many B-meson decays have been considered for observing effects originat- 
ing from physics beyond the Standard Model (SM). In general, the following 
classes of B decays are most sensitive to new physics effects: 1) Ab=l pro- 
cesses through penguin diagrams, 2) Ab=2 processes through box diagrams, 
and 3) tree-level processes mediated by exchange of a new particle. 

Processes such as b — > 57, belonging to the class 1), have been analysed 
but theoretical uncertainties hamper the observation of new physics 
signatures 0. Similar channels such as b —* sqq || and b — * s££ 0] also 
suffer from large theoretical uncertainties. Some other processes such as 
B — > t, representing the class 3), have been shown to be rather insensitive 
to a large class of new physics models ||. 

Rare decays, representing the class 2), can probe efficiently new physics 
effects, since for these decays the SM typically predicts extremely tiny 
branching ratios. The process b — ► ssd is a decay which is strongly sup- 
pressed in the SM. This decay can be produced in the SM by box diagrams 
(see Fig. |l]a) with an estimated branching ratio at a level lower than 10 -11 . 
The Minimal Super symmetric Standard Model (MSSM) introduces squark- 
gaugino (or higgsino) box diagrams (see Fig. fjjb), increasing the estimated 
branching ratio to 10~ 7 — 10~ 8 |J. The decay has also been studied in 
the context of several Two Higgs Doublet Models (THDM) {?]]. The studies 
have shown that in these cases, the branching ratio could be as high as 10 -7 . 
Supersymmetry with broken R parity provides another model with a signifi- 
cant enhancement of this decay (see Fig. ||c) || . These decays are tree- level 
processes (class 3) in our classification), and therefore the branching ratio 
could be even significantly higher than those predicted for the box-diagram 
processes. 
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a) Standard Model 



b) MSSM 



c) MSSM with broken R parity 



Figure 1: Feynman diagram contributions to the decay b — > ssd in various 
models: a) The Standard Model, b) MSSM, c) MSSM with broken R parity. 
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Three-body decays of a charged B, such as B — ► K K ir^ (either di- 
rectly or through a K*-resonance), were suggested as a clear signal for the 
process b — > ssd in Ref. ||. Recently, an upper limit of 8.79 • 10" 5 was set 
by the OPAL collaboration for the branching ratio BR(B ± — > K ± K ± 7r =F ) at 
90% confidence level @. 

In this paper, the direct decay B + — > K + K + 7r~ was considered [] in order 
to test the feasibility of observing these decays in the ATLAS experiment at 
the Large Hadron Collider (LHC). Similar analysis could be performed with 
a more general class of final states B 1 * 1 — ► K^K 1 * 1 + (no strange), including 
also K*-resonances, to increase the statistics. In sections 2-4 the analysis 
procedure and simulation results for ATLAS are described. In section 5 
the reach of other experiments is estimated, and section 6 summarises the 
paper. 

2 Event simulation 

The B + — > K + K + 7r _ decay was implemented in the Monte Carlo program 
PYTHIA 5.7 H in order to generate the signal sample. In the event gener- 
ation, b-quark pairs were produced in pp-collisions at y/s = 14 TeV either 
directly via the lowest order process, or via gluon splitting or flavour exci- 
tation. 

Events containing a B + meson were selected, and then the B + was forced 
to decay into a K + K + 7r~ final state. The associated b was forced to decay 
semileptonically into ^X, in order to satisfy the ATLAS level-1 trigger re- 
quirements for B hadrons (a muon with a p T > 6 GeV and \rj\ < 2.4)Q. The 
muon was not needed in the subsequent analysis as such. 

The ATLAS second level trigger for this hadronic B decay could be 
envisaged to be the presence of three charged particles with pr > 1.5 GeV, 
forming an invariant mass close to the B-meson mass. The detailed trigger 
rates have not been studied. 

For this feasibility study, a fast simulation program was used instead of 
a full GEANT simulation. The parametrisation was established by study- 
ing in detail the resolutions of the five helix parameters of the tracks in 
fully-simulated samples, including tails ( |]l0| , [jll|] ) . The smeared five helix 
parameters of the track and the corresponding covariance matrix were ob- 
tained, and a look-up table as a function of pt and rj was produced. The 
parametrisation was then applied to the four-momenta of the generated par- 
ticles. In case of pions, the parametrisation included a dependence on the 
decay radius as well, in order to be able to describe pions coming from the 
decay of long lifetime particles such as Kg. 

1 Charge conjugated states are implicitly included. 

2 Throughout this paper, the symbol px is used for the transverse momentum with 
respect to the beam direction, and i) for the pseudorapidity. 
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3 HADRON IDENTIFICATION 



3 Hadron identification 

The possibility for separating kaons, protons and pions enhances the observa- 



tion potential of many B-hadron final states in the ATLAS experiment 12]. 
Monte Carlo studies of K/tt separation using the signal shape information 
from the ATLAS Transition Radiation Tracker (TRT) have been previously 
reported in 

Recent test-beam data and detailed Monte Carlo simulations show that 
using the time-over-threshold information in the TRT data allows for an 
improved hadron identification, assuming that the TRT read-out would also 
provide the time of the trailing edge at low luminosity while preserving 
the output bandwidth requirements. The time-over-threshold method is 
described in detail in [14|. 



The time-over-threshold (ToT) for a single straw is defined as the width 
of the signal at the output of the low-threshold discriminator in the TRT 
front-end electronics. The ToT provides partial information on the particle 
energy loss in the straw gas, assuming that its dependence on the distance 
of closest approach of the track to the straw anode has been taken into 
account. 

The energy loss estimator ((At t)) is built on the basis of the individual 
ToT for all the straw hits on a given particle track, according to the proce- 
dure described in [O], In the TRT, on average, 35 straws will be crossed by 
particle tracks with > 0.5 GeV and \rj\ < 2.5. 

The expected K/tt separation as a function of the transverse momentum 
is shown in Fig. [2] in units of standard deviation. Without including any 
pile-up effects, the K/tt separation is predicted to be above one standard 
deviation for transverse momenta between 2 and 5 GeV, averaged over the 
full rapidity coverage (solid line), and above one standard deviation over a 
broader pj-range between 2 and 15 GeV at \rj\ = 0.3 (dotted line). 

In order to study physics processes in the ATLAS experiment, the mean 
and the sigma of the (At t) distributions, well described by gaussians, were 
parameterised for pions, kaons and protons as a function of px and r\ over 
the full acceptance region of the ATLAS TRT [15]. 
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Figure 2: Expected K/n separation in units of standard deviation in the 
ATLAS TRT as a function of transverse momentum (no pile-up effects in- 
cluded). The separation is shown as an average over the full rapidity cover- 
age (solid line) and at \rj\ = 0.3 (dotted line). 
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4 ANALYSIS 



4 Analysis 



4.1 Event selection 

The event selection cuts are summarised in Table [l]. The cuts on the trans- 
verse momentum of the particles and the loose cut on the B-candidate mass 
emulated the suggested second level trigger requirements. The other selec- 
tion criteria were based on the quality of the B-vertex fit, on the long lifetime 
of the B meson (the reconstructed B vertex was required to be separated 
from the primary vertex with at least 100 /jm in the transverse plane), on 
the rejection of events in which two of the particle pairs formed masses close 
to light resonance masses, and on the probability that the three particles 
formed a KK-zr combination. The last selection criterium is explained in more 



detail in Section \j.2[ The resolution of the decay length in the transverse 
plane was 72 /im. 



No. 


Selection 


Signal 


Background 






efficiency 


efficiency 


1. 


p T (tracks) > 1.5 GeV 


50.9% 


68.8% 


2. 


4 GeV < M(B) <6 GeV 


97.7% 


53.1% 


3. 


X 2 (triplet vertex fit) < 2 


83.6% 


68.4% 


4. 


p T (B) > 10 GeV 


80.1% 


22.6% 


5. 


Vertex detachment > 0.1mm 


58.3% 


0.3% 


6. 


V{dE/dx) > 0.1 


87.8% 


76.3% 


7. 


mj 3 and m| 3 > 2.5 GeV 2 


74.4% 


35.9% 


8. 


5.16 GeV < M(B) < 5.45 GeV 


91.2% 


5.8% 




Overall 


11.6% 


2.8 • 10~ 4 % 



Table 1: The signal and background efficiency of each selection cut. The 
cuts 1-8 were applied in sequence, and the efficiency of cut N is given relative 
to the remaining sample (N — 1 cuts applied). 

The overall signal efficiency was found to be 11.6%, while the background 
efficiency was 2.8-10 -4 %. The study of the background rejection was limited 
by the statistics of the simulated background sample (one million inclusive 
bb — > /i6X events, where //6 denotes the level- 1 trigger requirements for the 
muon). The generated background consisted of the default decay channels 
in PYTHIA 5.7. 



4.2 Use of the dE/dx information 



The selection criteria discussed in Section 4J. reduce the background by six 
orders of magnitude, while preserving about 12% of the signal, as it is shown 



4.3 Results 
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in Table [I]. In this section, the use of the dE/dx information is explained in 
more detail. 

For a charged particle track with a given and 77, the dE / dx ac tuai (j>T, f]) 
was simulated using the sigma and the mean provided by the parametrisation 
(see Section [|): 



dE 

dx 



(pr,v) = mean((A ToT ))(p T ,r/) + RND * <r((A ToT ))(p T , 77), 

actual 



where RND is a gaussian-distributed pseudo-random number. 

For any given triplet of particles, two positively and one negatively 
charged — candidates for the decay products of the B + - the % 2 dis- 
tribution was constructed according to : 



X 



dE I _ dEi I 
dx lexp dx I act 



where the index i labels the individual particles in the triplet and ^§|exp 
is the mean value of the dE/dx distribution for pions (if the particle had 
a negative charge) or for kaons (if the particle had a positive charge). The 
obtained x 2 probability distribution for three degrees of freedom is shown 
in Fig. ||. The background misidentification probability as a function of the 
signal efficiency is shown in Fig. when the cut on the dE /dx x 2 probability 
was varied. 

4.3 Results 

The number of signal events, passing the ATLAS level- 1 trigger, was esti- 
mated as: 



N *£L = ^(PP - bS - M6X) • 0r(b - B+) • Bt(B+ K+K+tO • J £dt, 
and the number of observed events as: 

/yobs _ Arprod 

JV signal — iv signal ' trec ' fc id> 

where the cross-section after the level- 1 trigger is o"(pp — > bb — > p6X) = 2.3/ib, 
Br(b — > B + ) = 39.7%, the integrated luminosity is J Cdt = SOfb^ 1 (corre- 
sponding to three years of LHC data-taking at the initial low luminosity), 
the signal reconstruction efficiency is e rec = 11.6%, the muon efficiency 
is e;d(^6) = 0.85, and the pion and kaon reconstruction efficiencies are 
e id (vr,K) = 0.90. 

The number of background events was estimated as: 
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4 ANALYSIS 





Signal efficiency 



Figure 3: The \ 2 probability 
distribution for three degrees of 
freedom for the signal (dotted 
line) and the background (solid 
line), after applying the proposed 
second level trigger requirements 
(see Sect. 2). Both distributions 
were normalised to unity. 



Figure 4: Background misiden- 
tification probability func- 
tion of the signal efficiency, when 
the cut on the dE/dx x 2 proba- 
bility was varied. The proposed 
second level trigger requirements 
were applied to the event sample 
first (see Sect. 2). 



Nl r g od = a(pp — bb — /x6X) • J Cdt. 

The number of produced background events was thus JVg° d = 6.9 • 10 10 for 
30 fb -1 . Using e rcc = 2.8 1(T 4 % and e id (/u6) = 0.85, e id (vr,K) = 0.90, the 
number of observed background events was N<* s = 1.2- 10 5 . 
An upper limit of 

Br(B + -> K+K+vr") = 3.4 • 10~ 7 

can be established at a 95% CL after three years of data-taking at the low 
luminosity. If one requires a signal significance of five standard deviations, 
a signal with a branching ratio of: 

Br(B + -» K + K+7T~) = 8.8 • 10~ 7 

can be observed with the same statistics. 

4.4 Limits on the R-parity violating couplings 

The b — > ssd decay rate induced by the R-parity violating couplings was 
estimated in Ref. @ to be : 
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r(b -»■ ssd) 



m b^QCD 
512(2vr) 3 m? 





3 


to 


3 


N5 


( 


/, -^32^21 


+ 


^12-^23 






i=l 




i=l 





where mi, is the b-quark mass, /qcd = (os(b)/a s (mp)) 24/ ' 23 , m p is the 
sneutrino mass and A' are dimensionless couplings. 

It was estimated in Ref. [|| that a fourth (or less) of the b — > ssd tran- 
sitions leads to B ± — > K ± K ± + (no strange) decay channels. Final states 
including K^K 1 * 1 + (no strange) can be produced via one or two excited 
kaons with respective proportions K*K* > K*K > KK. It was assumed, pes- 
simistically, that direct K ± K ± decay represents a third of the total K^K 1 * 1 
+ (no strange) decays, thus 1/12 of the b — > ssd transitions. Using this esti- 
mation with the 95 % CL ATLAS bound and, as in Ref. ||, = 4.5 GeV, 
/qcd = 2, mp = 100 GeV, t b + = 1.65 ps, a limit on the couplings will be: 





3 


2 


3 


2 




/] K32K2I 


+ 


K12K23 


< 5.3 • 10~ 5 




i=l 




i=l 





This limit should be considered as a rough estimate, obtained using 
very pessimistic assumptions on the relative decay probabilities in order to 
maximize the relative fraction of the decay mode B + — ► K + K + 7r~. The limit 
is nevertheless an order of magnitude better than the corresponding OPAL 
limit in Ref. ||]. Complementary measurements come from neutrino data 
[16] which give values for individual A' couplings for different neutrino mass 
scenarios. 



5 Comparison to the other experiments 

If the branching ratio of the decay B + — > K + K + 7r _ is in the range of the 
MSSM or THDM predictions (0(1O~ 7 )), the event yield of the PEP-II and 
KEKB B-factories is only a few events with an integrated luminosity of 30 
fb _1 to 100 fb -1 . Therefore these investigations would not seem feasible for 
BaBar and Belle in the R-parity conserving scenarios. 

The hadron colliders Tevatron and LHC have much larger bb cross sec- 
tions, which opens up the opportunity to study b — ► ssd transitions. For 
CDF it was assumed that the trigger efficiency is a ten times better than 
the ATLAS trigger efficiency, due to the possibility of triggering on purely 
hadronic final states. It was assumed that CDF has the same signal recon- 
struction efficiency as ATLAS. The K/ir separation capability is similar in 
the two experiments for transverse momenta above 1.3 GeV. ATLAS has a 
larger pseudorapidity acceptance than CDF, but on the other hand the final 
states are more central with less initial state gluon radiation at the smaller 
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center-of-mass energy. The estimated overall efficiency was consistent with 
the results on experimentally similar decays fi~7|| . Using these assumptions, 
one can estimate that CDF could set an upper limit of 7.9 • 10 -7 at 95% 
CL with 2 fb _1 . A signal significance of five standard deviations could be 
achieved if the branching ratio were 2.0 • 10 -6 . 

The LHCb experiment has the advantage of being able to trigger on 
purely hadronic final states, and having a superior K/ir separation thanks 
to its RICH detectors. Based on the results presented in Table 15.11 in 
Ref . Jl^l , it was estimated that LHCb could observe a five-standard-deviation 
signal if the branching ratio were 1.0- 10~ 7 , given the statistics of three years' 
running at the nominal LHCb luminosity of 2 • 10 -32 cm -2 s -1 . The 95% 
CL upper limit for the branching ratio would be 4.0 • 10 -8 . 

These results should be taken as crude estimations, which were based 
on the publicly available information on the detector and accelerator per- 
formance. 



6 Summary and outlook 

A feasibility study of reconstructing the decay B + — ► K + K + 7r~ in the AT- 
LAS experiment at LHC has been presented. The obtained 95% CL upper 
limit of 

Br(B + -> K+K+vr") = 3.4 • 10~ 7 

approaches the range of branching fractions predicted by MSSM or THDM 
scenarios. In R-parity violating models, branching ratios could be as large 
as 10 -4 . ATLAS could thus contribute in the measurements of some of the 
R-parity violating couplings. Given the upper limit above, the following 
limit can be set on the couplings: 
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2 


3 


2 




^i32^i21 


+ 


^12^23 


< 5.3 • 10" 5 








i=l 





The presented analysis of the B + — > K + K + 7r~ decay shows that the AT- 
LAS experiment will be able to set a new upper limit on the branching ratio, 
which will be more than two orders of magnitude lower than the present es- 
timate. The limit on the relation that constrains the A' couplings of the 
MSSM with R-parity violation will also be improved by an order of magni- 
tude. This analysis considered only the direct decay B + — > K + K + 7r~, but 
ATLAS should be able to increase the statistics by searching for final states 
with K*K* and K*K. 

Combining all the results, the LHC experiments will contribute signif- 
icantly to the search and measurements of physics beyond the Standard 
Model using the B + — > K + K + 7r~ decay channel. 
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